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Abstract 

Applying semi-classical Quantum Mechanics, the vacuum fluctua¬ 
tions within a star are determined, assuming a constant mass density 
and applying a monopole approximation. It is found that the density 
for the vacuum fluctuations does not only depend linearly on the mass 
density, as assumed in a former publication, where neutron stars up 
to 6 solar masses were obtained. This is used to propose a simple 
model on the dependence of the dark energy to the mass density, as a 
function of the radial distance r. It is shown that stars with up to 200 
solar masses can, in principle, be obtained. Though, we use a simple 
model, it shows that in the presence of vacuum fluctuations stars with 
large masses can be stabilized and probably stars up to any mass can 
exist, which usually are identified as black holes. 

1 Introduction 

The theory of General Relativity (GR) [1] has passed up to now all observa¬ 
tional tests. However, all these tests are for weak gravitational helds [2] (to 
which we also count the gravitational held in the Hulse-Taylor pulsar [3]), 
compared to the held strength near the Schwarzschild radius. 
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The standard GR predicts singularities, due to the black-hole solution. 
This solution has a singularity in the center and a coordinate singularity 
called the event horizon. No information below this event horizon can reach 
an external observer, even nearby. 

In our philosophical understanding, no theory should have a singularity, 
not even a coordinate singularity of the type discussed above. The appear¬ 
ance of a singularity hints to the incompleteness of a theory. While this is 
generally accepted for the singularity at the center, the discussion is going 
on whether this also applies for the coordinate singularity. 

In recent publications [4, 5] an algebraic extension of GR to pseudo¬ 
complex (pc) coordinates was proposed, denoted as pseudo-complex General 
Relativity (pc-GR). The pc-coordinates have the structure + ly^, 

with P = 1. Due to the property of = 1, variables with = dny^ have 
no inverse, thus the pc-variables form not a held but a ring. The absolute 
value squared | \ = X^*X^^ of these variables, with = x^ — ly^, are 

zero, though the variables themselves are different from zero. A modihed 
variational principle was proposed, which states that the variation of the 
action, which is now pseudo-complex, has to be proportional to a function 
whose absolute value is zero, called also a generalized zero. As a consequence, 
the Einstein equations change and, in the limit of y^ = 0, have the form 

, ( 1 ) 

where the the gravitational constant and the velocity of light were set to 1. 
The is the energy momentum tensor of the dark energy and the index 
A denotes this association. The same equation is obtained, adding in the 
standard GR by hand an energy-momentum tensor on the right hand side of 
the Einstein equations. 

Using semi-classical Quantum Mechanics (QM), where the background 
metric is hxed, it was shown that a mass changes the vacuum structure in 
its vicinity. The fluctuations increase toward the Schwarzschild radius and 
become inhnite there [6]. 

Because a quantized theory of gravitation does not exist yet, one relies on 
semi-classical QM, which does not permit to include the back-reaction of the 
vacuum fluctuation in a consistent way and is not applicable for strong grav¬ 
itational helds. Therefore, we proceeded to assume the presence of vacuum 
fluctuations, proposing a radial dependence such that for large distances. 
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compared to the Schwarzschild radius, the corrections to standard GR can 
not be measured yet. For the dark energy a fluid model was chosen, leading 
to corrections of the metric in the order of 1/r^. (In [7] other radial depen¬ 
dencies were investigated with similar qualitative results). This did permit to 
include the back-reaction to the metric and to solve the Einstein equations. 

The procedure reflects a general principle, namely: Mass not only curves 
the space but also modifies the vacuum structure, such that the metric itself 
changes. Because we used a fluid model, an additional parameter B = bwfi 
was introduced, where m is the mass in units of length and b is the actual 
parameter, measuring the coupling between the central mass to the vac¬ 
uum fluctuations outside the star. By an adequate choice of this parameter, 
namely that goo satisfies (l — ^ -f > 0, the b obtains a lower bound, 
such that there is no event horizonl 

In [7] the consequences for the case of the Kerr-solution were investigated. 
It was shown that the orbital frequency of particles in a circular orbit exhibit 
a maximum in the frequency which from there on is decreasing again toward 
lower radial distances. 

Observing Quasi Periodic Objects (QPO) in the accretion disk of galactic 
black holes [8, 9, 10, 11], assuming that they correspond to local bright spots 
in the disk following its orbital motion, one can deduce a radial distance using 
the prediction of the theory (GR or pc-GR). Measuring at the same time the 
redshift, using K — a lines, also has as a result a given radial distance. For 
a consistent theory, in both measurements the deduced radial distances have 
to be the same. And here is the problem: They agree in pc-GR but not in 
GR [12]. There is still a discussion going on if the QPO follows the orbital 
motion of the disk or are due to oscillations within the disk [13] provoked by 
the stellar partner, which would change the conclusions. Our argument that 
the QPO’s in accretions disks around galactic black holes correspond to local 
bright spots, following the motion of the disk, is based on their observation 
in the accretion disks around large masses in the center of galaxies, with no 
stellar partner nearby, and that the physics in both should be very similar. 

Another prediction of pc-GR is that the accretion disk appears brighter 
and exhibits a dark ring [14]. A model of a thin, infinitely extended disk [15] 
was applied. This dark ring is due to the maximum of the orbital frequency of 
particles in a circular orbit. At the maximum shear forces between particles 
in two neighboring orbits is small and thus little heat is additionally created. 

As one notes, pc-GR provides definite observable predictions, which will 
be observed in near future. 
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A further application of pc-GR is given in [16], where the question was 
investigated if due to the presence of dark energy, neutron stars with larger 
masses than 2-3 solar masses are possible. Due to the lack of knowledge on 
how the dark energy couples to the mass distribution within a star, a simple 
ansatz was used, namely that the dark energy density ev is proportional to 
the mass density e^- The proposed coupling reads ca = 06^, with a < 0. As 
a result neutron stars up to 6 solar masses were obtained and their stability 
was proven. Larger masses could not be realized, because the coupling turned 
out to be too strong near the surface. We attributed it to the proportionality 
between the two densities. 

In order to get a better estimate of the coupling between the dark energy 
and the mass distribution, we are lead to consider semi-classical QM with 
mass present. It refers to the interior of normal stars and known neutron 
stars. The suggested coupling will be treated as a phenomenological model. 

The outline of the paper is as follows. In Section 2 we present the main 
points and formulas in the semi-classical treatment, hrst outside a star and 
then within a star, using the monopole approximation which permits a sim¬ 
pler treatment. The Schwarzschild metric is considered as the background, 
i.e., we discuss non-rotating stars. Two different kinds of mass distributions 
are considered. The hrst is a constant mass density within the star and the 
second one is a density which falls off such that it simulates realistic density 
distributions, except near the surface of a star. We will show that the cou¬ 
pling between the dark energy and mass density in the interior region has 
to diminish toward the radius of the star, even for a constant mass density. 
The models discussed are simplistic but reveal the key properties. 

In the same section, a simple radial dependence of the dark energy to the 
mass density is proposed and used in Section 3 to calculate the masses for the 
stars. We show that now up to 200 solar masses are possible. Higher masses 
could not be obtained due to numerical difficulties. The result indicates 
that it is a matter of the correct coupling that ’’neutron stars” of arbitrary 
mass can be created, or in different words, that the large masses observed in 
the center of nearly every galaxy are ’’neutron stars”, which are very black 
and rather simulate a black hole. The main motivation is to give a prove 
of principle that stars can acquire arbitrarily large masses, without forming 
a black hole. (Though, it is doubtful that such objects still resemble the 
neutron stars as we know them.) 

In Section 4 conclusions are drawn. 
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2 Semi-classical Quantum Mechanics and Vac¬ 
uum fluctuations 

A very good introduction to semi-classical Quantum Mechanics is given in 
[17]. The main goal is to determine the vacuum fluctuations in a flxed back¬ 
ground metric, for example a Schwarzschild metric. The physical quantity of 
interest is the expectation value of the energy-momentum tensor T^j,, which 
is a local quantity and is preferred versus the observation of particle numbers. 
The measure of particle numbers is observer dependent, while transforms 
under a relativistic transformation, relating one system to another equivalent 
one in a well defined way. The energy-momentum tensor has to be regularized 
and renormalized, where the methods are explicitly presented in [17]. 

The regularization/renormalization of the energy-momentum tensor is 
still quite involved and approximate methods have to be applied in order to 
determine its expectation value in 4-space-time [6, 18]. A simpler approach 
is followed invoking the monopole approximation [17, 19]. The monopole 
approximation consists of assuming a spherical symmetry and restricting the 
length element squared to the time and radial component, i.e. 


ds^ = e'' dt^ — e^dr^ 


( 2 ) 


with goo = and gu = e^. Defining the tortoise coordinate [19] via 


_ e^-‘'dr^ 


(3) 


leads to 


= C{r) {dt^ — d^"^) , (4) 

with C{r) = = goo- Because we study a time-independent spherical metric, 

the factor C depends on the radial distance only. 

Defining further 


f ^ , X = t ^ ^ t i (x+ + X ) , ^ = ^ (a:+ - X ) ,(5) 
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leads to the line element 


= C{r)dx~^dx 


( 6 ) 


In this 2-dimensional space the expectation value of can be deter¬ 
mined readily [17, 19]: 
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Transforming the coordinates x± to t and r, we obtain the alternative 
expression 
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Using the Tolman-Oppenheimer-Volkov (TOV) equation [20], which relates 
Z 2 , A and their derivatives in r to the density, pressure and the derivative of 
the pressure, we obtain 
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The 4-dimensional result is approximated by 


(9) 



( 10 ) 
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with = T^y. This is the monopole approximation. The advantage of it 
is its simplicity, but deviations from the exact result are to be expected. 

For the Schwarzschild solution, in [19] the expectation value of the energy- 
momentum tensor components in the exterior region of a star was determined 
in two dimensions. After the transformation to the coordinates t and r, one 
obtains for the 4-dimensional expectation values outside the mass distribu¬ 
tion, defining = 2m (the Schwarzschild radius) and skipping the index 


4D”, 



( 11 ) 


where we have transformed to the mixed components of the tensor (one 
component below and the other above), which permit a direct comparison to 
the density, e.g., T/ = e\. The index A refers to its interpretation as dark 
energy. The result agrees in structure with [6], though the factors are not 
the same. 

As can be seen, the fall-off of the density is in leading order proportional 
to 1/r®. This was also obtained in [6] but with a different factor, also in the 
correction terms. However, as can be noted in (11), a singularity appears 
at the Schwarzschild radius. There the vacuum fluctuations tend to infin¬ 
ity, rendering the semi-classical approach invalid. The qualitative conclusion 
drawn from this is an increase toward the central mass with some 1/r” depen¬ 
dence. The origin of this coordinate singularity at the Schwarzschild radius 
is to neglect the back-reaction on the metric. 

This observation leads to the proposal of the pseudo-complex General 
Relativity (pc-GR) [4, 5]. In this modihed theory of General Relativity, an 
additional term on the right hand side of the Einstein equations is required, 
which we associate with the energy-momentum tensor of the vacuum fluctua¬ 
tions. For the energy-momentum tensor a model of an ideal asymmetric fluid 
is applied. The advantage of this is that the back-reaction of the vacuum fluc¬ 
tuations on the metric can be determined. However, due to its phenomeno¬ 
logical nature, an additional parameter appears. In its first version [5] the 

density of the dark energy falls off as 1/r^ and leads to goo = — ^ + j • 

Another interest is to investigate the vacuum fluctuations within a star, 
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which requires the knowledge of the coupling of the dark energy to the mass 
distribution. In [16] a simple coupling model was applied in which the ratio 
of the dark energy density to the mass density is a constant. In what follows 
we will estimate, using simple assumptions, the vacuum fluctuations as a 
function in the radial distance, which will lead to a better proposal for the 
coupling, finally applied in the next section. 

For the ansatz of the mass distribution two models are considered: i) A 
constant mass distribution within the star, which is the easiest to treat, and 
ii) a density which simulates a realistic behavior. These distributions are im¬ 
plemented by hand and not derived from the Tolman-Oppenheimer-Volkoff 
(TOV) equations, required by a consistent approach. We are, however, not 
interested in an exact description but rather in obtaining an idea and moti¬ 
vation for a phenomenological ansatz of the dark energy density as a function 
in r. This will then be applied in the next section. 

2.1 A constant mass distribution 

The metric for a constant mass distribution Sm = is derived in [20], i.e.. 



-r^ - sm\d)dip^) , (12) 

with R as the radius of the star and f'l = 

With (8) and (9), using the monopole approximation (10), we obtain for 
the expectation value of the T/ and component 


£a 


1 



1 (gp -F 3prf 
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Pr = = 7 ^ (^0 + 

r2(£o-Pr-)(£o + 3p^)\ 

+ 2 -— —) 

The prime in p], of the matter pressure refers to a derivative in r. We do 
not present the part for the matter density and the derivative of its pressure, 
because we are interested in the dark energy density part only. 













Suppose that inside the star the fluid is isotropic. Then Pr = p and using 
the TOV equations (see [16]) for the case of a constant mass density, we 
obtain [20] 



From now on, only e\ will be considered. Substituting (14) into (13) 
leads for the dark energy density 



The factor in the numerator was rewritten as follows: With m = ^SqR^, 
resolving for R^, we obtain with ^0 = ^ (see the definition below (12)) the 
value ^0 = g^- In fhe factor in (15), this is substituted only for one Sq. 
In this manner we show the explicit dependence of Sa to eo- 

The density, as given in (15) is expanded up to leading order in r^. In 
Fig. 1 the qnadratic approximation is compared to the exact relation (15). 
Note that the major differences only occur near the surface. However, the 
approximation gets worse for lower radii R. 

In conclusion, the dark energy density £a is not only proportional to the 
constant mass density with a constant factor, but there is an additional r- 
dependence, which can be approximated for R not too small by a factor of 
the type (1 — (3r^). The constant (3 depends in (15) on the constant energy 
density and on the radius of the star. Assuming that the same relation 
of the coupling is valid for a matter density which also does depend on r, 
the coupling according to (15) should acquire an even more complicated r- 
dependence. However, in order to keep the matter simple, the ansatz which 
depends on r^, should suffice for a start. 

To resume, the considerations in this section suggests the following ansatz 
for the coupling of the dark energy to the matter density: 


£a = 



(16) 
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Figure 1: The figure compares the radial dependance of the dark energy 
density for a star with constant matter density (see eq (15)), radius R = 
ym 4.6m and Eq = 5767rrQ with its second-order approximation. 
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where a is a factor as used in [16] and rg is an additional parameter, not to 
he confused with the parameter used in the model of a constant mass density. 

2.2 A non-constant mass distribution 

Let us assume that the mass density varies as 




m 


^0 


i-t, 



(17) 


which simulates realistic calculations [16]. As a length element squared we 
use 


= e'^dt‘^ — e^dr"^ — {dd^ + sin^'ddyj^) , (18) 

where the metric for the interior of the star is proposed to behave as 



In order that this connects smoothly to the Schwarzschild metric outside the 
star, the following condition has to be fulfilled: 



2m 
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( 20 ) 


Following the same steps as described further above, we 
dark energy density at 


arrive for the 
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The obtained r-dependence of the dark energy density on the matter 
density is now, as expected, more complicated. Nevertheless, the main (1 — 
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6r^)-dependence is reproduced, with the correction factor proportional to 
Only near the surface (r approximately R), this correction is large and 
may be attributed to the ansatz of the matter density, i.e., it is expected to 
be not very good there 

Again, we arrive at the structure similar to the ansatz (16). Thus, in the 
next section we use this ansatz to investigate if stars with larger masses than 
those obtained in [16]. 

3 Large mass stars 

The theory of neutron stars within the pseudo-complex General Relativity is 
published in [16], assuming a purely linear coupling between the dark energy 
density and the corresponding quantity for the baryonic counterpart Em- 

£\ = OcEm ( 22 ) 

The coupling parameter was chosen to lie in the range (—1,0) and under 
such conditions, stable solutions were obtained having masses for the bary¬ 
onic component of up to [6 — 7]Mq. In the present work, we repeat such 
calculations using instead of (22) the new coupling relation (16). 

As can be appreciated from Fig. 2, the new coupling relation allows 
the study of solutions with values of a < —1 whose masses can acquire 
considerable high values of up to 200 solar masses, in comparison with the 
previous work. This can be understood as a consequence of the stronger 
fall-off of the new relation toward the radius of the star that causes the dark 
energy repulsion near the surface to be weaker than the pure linear model 
(22), while it remains strong near the center. Thus, larger stellar masses can 
be maintained without collapsing and the repulsion is not strong enough to 
evaporate the surface. 

The presence of these more massive solutions can be appreciated in Fig. 
3 where already one can identify those which have ^ < 0 and therefore 
are unstable against small perturbations. Unfortunately, the stability of the 
solutions with ^ > 0 cannot be completely assured by simple criteria as 
those employed in [16] and a full perturbative method needs to be carried 
out. Numerical instabilities inhibited a more profound study. Checking the 
previous, necessary condition however tells us that at least such massive 
objects could in principle exist. 
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Figure 2: Mass of the baryonic component vs. radius for different values 
of the coefficient a. The central A-term pressure p\c has been hxed to 
Pemc=io-'^eoj fh® tq parameter to 100 km and the central baryonic energy 
density Smc between [4 • 10“^ — 10]eo- 
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Figure 3: Mass of the baryonic component vs. its central energy density for 
different values of the coefficient a. The central A-term pressure p\c has been 
hxed to Pe^^=io- 2 £o, the ro parameter to 100 km and the central baryonic 
energy density Emc between [4 • 10“^ — lOj^o- Stability can be proven, as 
explained in [16], investigating how the curve bends for large central density. 
However, we could not proceed because data on the baryonic properties were 
not available for that region. 
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Figure 4: Mass of the baryonic component vs. radius for different values 
of the parameter tq. The central A-term pressure pac has been hxed to 
Pemc=io-'^£o^ coupling parameter a to —1.2 and the central baryonic energy 
density Smc between [4 • 10“^ — 10]eo- 
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Within the current study, an additional parameter tq is introduced. As 
can be seen in Figs. 4 and 5, larger values of vq translate into larger values of 
the mass corresponding to the baryonic component. The value of ro cannot 
be a priori determined and therefore four values covering the whole range in 
which the radii of the solutions lie^ have been chosen. 



Figure 5: Mass of the baryonic component vs. its central energy density for 
different values of the parameter tq. The central A-term pressure pac has 
been fixed to Pe^^=io-2eo’ coupling parameter a to —1.2 and the central 
baryonic energy density Smc between [4 • 10“^ — lOj^o- 


^Please note that the value chosen for ro will modify the stellar radius. 
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4 Conclusions 


Using semi-classical Quantum Mechanics, the mass-induced vacuum fluc¬ 
tuations outside and inside a mass distribution were determined within a 
monopole approximation. Outside a mass distribution the dominant part of 
the dark energy density falls off as which is in qualitative agreement to 
calculations in four dimensions [6]. This finding justifies the ansatz for the 
dark energy density, as used in the pc-GR [4, 5]. 

Inside the mass distribution similar calculations lead to a deviation of 
the linear coupling, as assumed in [16]. The encountered dependence of the 
coupling between the dark energy density and the mass density includes an 
extra factor, depending on the radial distance r. This finding permits a new, 
phenomenological ansatz for the coupling, which allows a large dark energy 
density in the center of a star, therefore sustaining a large mass, and a low 
density near the surface, thus not evaporating the outer layers of the star, as 
found in [16]. 

The discussion presented is a prove of principle that arbitrary large masses 
of stars are possible, which indicates that the so called black holes may be 
in fact large stars, though due to the strong gravitational held the physics at 
the surface and in its interior is expected to change significantly, compared 
to known neutron stars. 
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